NEW COMPUTABLE NECESSARY CONDITIONS FOR 
THE REGULARITY THEORY OF OPTIMAL 
TRANSPORTATION 

PAUL W.Y. LEE 

Abstract. We give new computable necessary conditions for a 
class of optimal transportation problems to have smooth solutions. 



1. Introduction 

Let fi and be two Borel probability measures on the manifold 
M and let c : M x M — y M be a cost function. In the optimal 
transportation problem, one looks for a Borel map which minimizes 
the following total cost among all Borel maps (p : M — > M which 
push /i forward to u: 

c(x, ip{x))dfi{x) 

Here the push forward ip^:fi of a measure yU by a Borel map is the 
measure defined by <^*/i(t/) = /i((^~^(f/)) for all Borel sets U C M. 

Under some mild assumptions on the cost c and the measures n, u, 
the above problem has a unique solution [U [T71 [3l O [H |8]. This 
unique solution is called the optimal map. There are various recent 
breakthroughs in understanding the regularity of the optimal map 
[m [IHl [131 HH HI]- The most important one is the introduction of 
a geometric object called the Ma-Trudinger-Wang (MTW) curvature. 
After the work of fTEl \TE\ IT3] . it is clear that certain non- negativity 
condition on the MTW curvature, called MTW condition, is necessary 
for the regularity theory of optimal maps. However, if the cost is not 
given by an explicit formula, then it is very hard to compute the MTW 
curvature and the MTW condition. When the cost c is given by square 
of a Riemannian distance, the following is the only known computable 
condition which is necessary for the MTW condition. 

Theorem 1.1. [IS] Let d be a Riemannian distance function on the 
manifold M and assume that the cost c is given by c = d^. Then the 
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MTW curvature satisfies 

MTW{u,0,w) = K{u,w){\u\'^\w\'^ - {u,wf) 

where K{u,w) is the sectional curvature of the plane spanned by u and 
w. 

In particular if the cost c satisfies the weak MTW condition (A3w) , 
then the sectional curvature is non-negative. 

The purpose of this paper is twofold. First, we consider cost func- 
tions arising from natural mechanical systems. More precisely, let (■, ■) 
be a Riemannian metric on the manifold M and let | ■ | be the cor- 
responding norm. Let V : M — > M be a smooth function on the 
manifold M, called the potential, and let L : TM — )■ M be the La- 
grangian defined by L{x,v) = — V{x). The cost functions that 
we are interested in are given by 

(1.1) c{x,y) = mf [ Lijit),jit))dt, 

Jo 

where the infimum is taken over all smooth curves 7(-) satisfying 7(0) = 
X and 7(1) = y. 

In the first part of the paper, we give computable necessary condi- 
tions for the cost defined in (11.11) to satisfy the weak MTW condition 
(A3w) (see Theorem 14. ip . The following is a simple corollary of The- 
orem |4]T1 

Theorem 1.2. Let x be a maximum point of the potential V. Assume 
that V also satisfy the following 

VK- = and RessV^ = 0. 

Let u and w be two tangent vectors based at x. Then the MTW curva- 
ture for the cost c defined in U.l\) satisfies 

MTW{u, 0, w) = {R{w, u)w, + ^ (V^V.VK, u) . 

In particular if the cost c satisfies the weak MTW condition (A3w) , 
then the following holds 

{R{w, u)w, u) + ^ (V^ V^VK, u) > 

for all orthogonal pairs [u, w) of tangent vectors {u, w) = 0. 

We remark that the condition Hessl4 = is not completely necessary 
and it can be replaced by a more complicated condition (see Theorem 
14. ip . Note that when the potential = 0, the cost c is given by the 
square of the corresponding Riemannian distance d and Theorem 11.21 
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reduces to Theorem ll.li As a corollary of Theorem 11.21 we have the 
following. 

Corollary 1.3. Let A be a n x n matrix satisfying 

{{Au, w) + {Au, w) f + 2 {Au, u) {Aw, w) > 

for a pair of vectors {u, w) in M" which are orthogonal {u, w) = 0. 
Let V : M" be a potential satisfying 

V{x) = - {Ax,xf + 0{\xf) as \x\ ^ 

and let L be the Lagrangian defined by L{x, v) = ^I'wp — V{x) where \ ■ \ 
is the Euclidean norm. 

Then the MTW curvature for the cost c defined in / fi.ij) does not 
satisfy the weak MTW condition (A3w) . 

In the second part of the paper, we focus on the case c = d'^, where 
c? is a Riemannian distance on the manifold M. We go beyond Theo- 
rem 11.11 and consider higher order necessary conditions for the MTW 
conditions. More precisely, according to Theorem II. H the Riemann- 
ian manifold M necessarily has non-negative sectional curvature if the 
cost c = d"^ satisfies the MTW conditions. However, when the sec- 
tional curvature is only non-negative. Theorem 11.11 does not tell us 
anything about the MTW conditions near where the sectional curva- 
ture vanishes. To understand the MTW conditions near these points, 
we consider the higher order Taylor expansion of the MTW curvature 
in the f-variable. If we assume that the sectional curvature K{u, w) 
of the plane spanned by u and w vanishes, then the zeroth order term 
in the Taylor expansion of the MTW curvature MTW{u,v,w) in v at 
the origin vanishes by Theorem 11.11 Therefore, if the MTW curvature 
satisfies MTW{u,v,w) > for all small enough v, then necessarily the 
first order term in the Taylor expansion vanishes and the second order 
term is non- negative. As a result, we get new necessary conditions for 
the cost to satisfy the weak MTW condition (Theorem I5.ip . When 
the manifold is two-dimensional, the conditions are simplified and give 
the following simple result. 

Theorem 1.4. Assume that M is a two dimensional Riemannian man- 
ifold with Riemannian distance function d. If the cost c = d^ satisfies 
the weak MTW condition (A3w) , then M has non-negative Gauss 
curvature and the Riemannian curvature R satisfies 

3 {{VuV^R){w, u)w, uf < 2 {{VlR){w, u)w, u) {{VlR){w, u)w, u) 

for each pair {u, w) of orthogonal vectors {u, w) = which spanned a 
plane with zero sectional curvature (i.e. {R{u,w)u,w) = 0). 
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As an example, we consider the two dimensional Euclidean space 
equipped with the metric 

(1.2) (u, t;) = e^^^^'^u ■ 

where u ■ w denotes the usual dot product and /(x, y) = x^y + ax'^y'^ + 
xy^ + O'iy'^- 

For these Riemannian metrics, the Gauss curvature is zero at the 
origin and nonzero everywhere else if a < —3. As a result of Theorem 
11.41 we get the following. 

Theorem 1.5. Assume that the Riemannian distance d is defined by 
the Riemannian metric (■, ■) given in U.2\) . If the cost c = d"^ satisfies 
the weak MTW condition (A3w) , then 



a<- 




2. Background: The MTW Curvature 

In this section, we will review some basic facts about the optimal 
transportation problem and the definition of the Ma-Trudinger-Wang 
(MTW) curvature. The assumptions in the theorems stated in this 
section are simplified to avoid heavy notation. The corresponding the- 
orems with relaxed assumptions can be found, for instance, in [19]. 

Let (■, ■) be a Riemannian metric on a manifold M and let \^ : M — >^ 
R be a smooth function which is bounded above. Let L : TM — )■ R be 
the Lagrangian defined by 

L{x,v) = 2l'^P ~ V{x). 
In this paper, we are mainly interested in the cost c defined by 

(2.1) c{x,y) = mf [ L{^{t),^{t))dt, 

Jo 

where the infimum is taken over all smooth curves 7(-) satisfying 7(0) = 
X and 7(1) = y. 

Curves t ^ '-/(t) which achieve the above infimum are called curves 
of least action and they satisfy the following equation (see [T2] ) 

(2.2) = -Vy,. 

Here we abuse notation and denote the covariant derivative by dt- The 
same convention will be used throughout this paper. 
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If t 7(t) is a curve of least action with initial velocity f , then the 
c-exponential map exp*^ is defined by 

exp'^iv) = 7(1). 

Note that, unlike the Riemannian case, t h- )■ exp'^{tv) is not a curve of 
least action in general. 

Let /i and u be two Borel probability measures with compact sup- 
ports on the manifold M. We recall that the optimal transportation 
problem is the following minimization problem: 

Find a Borel map which minimizes the following total cost among 
all Borel maps cp : M — y M which push fi forward to u: 

c{x, ip{x))dfi{x) 

M 

Here the push forward of a measure /i by a Borel map ip is the 
measure defined by (p^fi{U) = n{(p-^{U)) for all Borel sets U C M. 

Theorem 2.1. Suppose that the cost c is given by 1(2. 1\) and the mea- 
sure fi is absolutely continuous with respect to the Lebesgue measure. 
Then there is a solution ip (called the optimal map) to the above optimal 
transportation problem which is unique fi-almost everywhere. More- 
over, there exists a Lipschitz function / : M — )■ M such that the unique 
optimal map if is given by 

ifix) = exp'^(V/(x)). 

Next, we discuss the main object of this paper, the Ma-Trudinger- 
Wang (MTW) curvature. Let u, v, and w be vectors based at the point 
X. The MTW curvature MTW is defined by 

3 

MTW(m, V, w) = — d'}d'!c(a(t), exp%v + sw)) 

2 s=t=o 

where a is any curve with initial velocity u (i.e. dfcr = u). 

t=o 

Finally, we can state the MTW conditions. Let O be the set of all 
pairs of points (x, y) contained in the product M x M such that 

(1) there exists a unique curve of least action 7 satisfying 7(0) = x 
and 7(1) = y, 

(2) the map dexp'^ is a submersion at 7(0). 

It is known that the cost function c is smooth on the set O (see, for 
instance, |T2|) and the MTW curvature is well-defined. Let O be the 



subset of all initial velocities dt'~f in the definition of O. Then the 

t=o 

weak MTW condition is given by the following: 
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The cost c satisfies the weak MTW condition (A3w) on a subset ^A 
of O if 

MTW{u,v,w) > 

on the set 

{{u,v,w)\v G O, {x,exp''{v)) E M,{u,w) = 0}, 

The relevance of these conditions to the regularity theory of optimal 
maps can be found in[T6lll3l|l2[IIl[l8l|15l[71|8l[i[T0l|6]. 



3. The Ma-Trudinger-Wang curvature and the 
riemannian curvature 

In this section, we give a formula for the MTW curvature in terms 
of the change in the Riemannian curvature and the Hessian of the 
potential along curves of least action. Before stating the precise result, 
let us introduce the following notations. Let u, v, and w be tangent 
vectors based at the point x and let r 7s (t) be the curve of least 

action with initial velocity v + sw (i.e. dr'^s = v + sw). Let r i— )■ 

T=0 

Us{t) be the parallel translation of the vector u along the curve r i— )■ 
7s (r). Let r i— )• Js(r) be a vector field defined along the curve r i— ^ 
7s (r), called Jacobi field. It is defined as the solution of the following 
Jacobi equation 

d^J + Ridri, J)dri + Hess\/^(J) = 0. 

We assume that the Jacobi field J(-) also satisfies the following bound- 
ary conditions Js(0) = m, Js(1) = 0, and Js{t) 7^ for all time r in the 
interval (0, 1). 

Theorem 3.1. The MTW curvature is given by 
MTW{u,v,w) 

(3.1) _ 3 r r 

~ 2 X Jo 

Proof. The Jacobi field J satisfies the following Jacobi equation 
d'^J + R{dri, J)dri + Hess1/^(J) = 0. 
It follows that 



dl {R{dr-f, J)dri + Hess\/^( J), U) drdf 



dr {drJ, U) + {R{dri, J)(9,7, U) + (HessV;(J), [/) = 0. 
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If we integrate with respect to the variable r, then the above equation 
becomes 



{U.drJ) 



T = 



{U.drJ) _+ / (/?(9,7, J)9,7 + Hess1/^(J),t/)c/r 



dr {U, J) 



{Ridr-f, J)dri + Hess1/^(J), U) dr. 



Now if we integrate again with respect to r and use the boundary 
conditions for J, then we have 



(3.2) 
By 



{U.drJ) 

r=0 

= - + 



1 rf 



^0 



(i?(9,7, J)dr-i + HessV:,(J), U) drdf. 



Theorem 3.1], we know that the MTW curvature is given by 



MTW(u, V, w) = -dl {U, drJ) 



T = S=0 



The result follows from this and fl3.2D. 



□ 



4. Zeroth Order Condition for Natural Mechanical 

Actions 

In this section, we give the proof of the following main theorem. 

Theorem 4.1. Let x he a maximum point of the potential V . Let u 
and w be two tangent vectors based at x. Then the MTW curvature for 
the cost c defined in ^2. 1\) satisfies 



MTW{u,Q,w) 



3 
2 



1 /-f 



Jo 



2 {R{drW, u)drW, u) 



+ ^HessK(M), / R{drW, w)udr ) + (V^V^VK, u) drdf 



s=0 



where u and w satisfies the following linear ordinary differential equa- 
tion 



dlu = - 

w satisfies the initial conditions 



-HessK-fw) 



w 



drW 



r=0 



T=0 
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and u satisfies the boundary conditions 

0, M^O z/0<r<l. 



u 



T=0 



T = l 



Let us first give the proof of Theorem 11.21 



Proof of Theorem \l.S[ Since HesSa;^ = 0, we have w = tw and u = 
(1 — t)u. If we substitute this back into the formula for MTW(m, 0, w) 
in Theorem I4.H then we have 

MTW{u,0,w) 

o rl j-f 

2{l -t){R{w,u)w,u) 



s=0 



2 ./o JO 



+ T\l-T){VlVuVV^,u)dTdf 



{R{w, u)w, m) + — {VlV^VV,, u) . 



□ 



Proof of Theorem \4-l\ First let us note that the cost c is smooth at the 
point {x^x). Indeed, since x is a maximum point of the potential V, 
the constant curve 7(-) = x is the unique minimizer satisfying 7(0) = x 
and 7(1) = X. Let r 1— )• J(r) be a vector field defined along 7 which 
satisfies the Jacobi equation 

dlJ + R{dri, J)dri + HessV;(J) = 0. 

Since 7 = 0; and 9^7 = 0, it follows that 

9^J + Hess\4(J) = 0. 

The point x is a maximum point of the potential V , so the Hessian of 
the potential HessV^ is non-positive definite. Therefore, if J satisfies 

the boundary conditions J 



J 



0, then J = 0. It follows 
has full rank at the 



from [12, Theorem 7.6] that the map ciexp 

origin. Therefore, the cost function c is smooth at the point (x, x) by 
[121 Theorem 7.7] (see also Section [2]). 

Let us first introduce some notations. Let r 1— )■ 7s,t(T) be the curve 
of least action with initial velocity tv + sw. Let r 1— )■ Us,t{T) be the 
parallel translation of the vector u along the curve r 1— t- 7s,t(T). Let 
r I—)- Js,t{j) be the Jacobi field defined along the curve r 1— i- 7s,t(T) 
which satisfies the conditions Js,t(0) = m, Js,t(l) = 0, and Js,t{j) 7^ 
for all r in the interval (0, 1). 
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If we expand the term dl {R{dr'j, J)dr'y, U) in f l3.ip . then we have 

dl{R{drl.J)drl.U) 
= {{dlR){drl, J)drl, U) + {R{dld,l, J)drl, U) 

+ {R{d^^,J)d^^,dlU)+2{{dsR){dsd^i,J)d^^,U) 
(4.1) + 2 {{dsR){d^i, dsJ)dri, U) + 2 ((9,/?)(9.7, J)9,9.7, t^) 
+ 2 ((9,i?)(9.7, J)drl, dsU) + 2 (i?(9,9.7, dj)d^^, U) 
+ 2 (/?(9,9,7, J)dsdrl, U) + 2 (i?(9,a,7, J)drl, d^U) 
+ 2 (i?(5.7, d,J)dsdri, U) + 2 (i?(5,7, dsJ)dri, d^U) 
+ 2 (/?(9,7, J)dsdr-f, dsU) . 
By (1) of Lemma [8. H (14.11) simphfies to 



dl{R{d^^,J)d^^,U) 



s=t=0 



2{R{d,dr-i,J)dAi,U) 



s=t=0 



If we apply (2) of Lemma 18. H (1) of Lemma 18.21 and Lemma 18. 3[ 
then the above becomes 



(4.2) 9^ (R(dr-f, J)dr-f, U) =2 (Rid^w, u)drW, u) . 

s=t=0 

If we expand the other term (9f (Hessl^(J), U) in (13. ip . then we have 
{liessV^{J),U) = (Hess\/^(J),9f?7) 

(4.3) +(92(Hessl^^(J)),f/) 

+ 2(a,(Hess\/^(J)),9,f/). 
By (1) and (2) of Lemma [8l2l (14.31) becomes 



92 (Hess\/^(J),f/) 



s=t=0 



92(Hess\/^(J)),f/) 
+ ^ IIess\4(-u), / R{drW,w)udT 



s=t=0 



By (1) of Lemma [8. 3[ it follows that 



(Hess\/^(J),t/) 



s=t=0 



+ ^IIessV;(M) 



R{drW, w)udT 



Finally we combine this with (14. 2 p and Theorem 13.11 to finish the 
proof. □ 
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5. Higher Order Conditions in the Riemannian Case 

In this section, we consider the first and the second order terms of 
the MTW curvature in the f-variable. More precisely, we will prove 
the following second main result of the paper. 

Theorem 5.1. Let d be a Riemannian distance function on the man- 
ifold M. Assume that the cost function c is given by c = d"^ Then the 
MTW curvature satisfies 



dtMTW(u,tv,w) 

t=o 

= - ((V^i?)(w, u)v, u) + - {{W^R){w, u)w, u) 



and 



d:MTW(u,tv,w) 

t=o 

= ((V^i?)(t;, u)v, u) - - {R{v, u)u, R{v, w)w) 

4 2 

+ — {R{v, u)v, R{w, u)w) + - ((Vi,V^-R)(w, u)v, u) 
15 5 

+ T7T ^)^' T (^(^' ^)"' ^(^' ^)^) 

10 5 

4 

H ((R(w, u)v, R(w, u)v) + (-R(f , u)w, R(w, u)v)) 

15 

+ - {{R{w, u)v, R{v, w)u) + (-R(f , u)w, R{v, w)u)) . 

The proof of Theorem 15. II will be postponed to Section[71 For the rest 
of this section, we will state and prove the consequences of Theorem 
O 

Theorem 5.2. Let d be a Riemannian distance of non-negative sec- 
tional curvature on the manifold M. Assume that the cost c = d"^ 
satisfies the condition (A3w) . Then for each pair [u, w) of orthogonal 
vectors {u, w) = for which the plane spanned by u and w has zero 
sectional curvature, we have 



(5.1) 



{{yu,R){w,u)v,u) = 
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and 

G{u,v,w) := ^ {{VlK){v,u)v,u) - ^ {R{v,u)u, R{v,w)w) 

2 1 

+ 3 {{V,V^R){w, u)v, u) + — {iVlR){w, u)w, u) 
5 10 ^ ' 

4 

H ((R(w, u)v, R(w, u)v) + (Riv, u)w, R(w, u)v)) 

15 

H — {{R{w, u)v, R{y, w)u) + {R{y, u)w, R{y, w)u)) > 0. 
3 

Proof. We extend the vectors u,v,w to vector fields U,V,W, respec- 
tively. Moreover, we assume that U, V, W are constant vector fields 
in a geodesic normal coordinate neighborhood of the point x. By as- 
sumption, the function {R{W, U)W, U) has a minimum at x. It follows 
that 

{{y^R){w,u)w,u) = V {R{W,U)W,U) = 0. 

This proves the first equality. For the inequality involving G, we 
need to show that R{u,w)u = R{w,u)w = 0. Indeed, we know that 
R{u, ■)u is a symmetric operator. Since the manifold has non-negative 
sectional curvature, R{u, ■)u is non-negative definite. We also have 
{R{u, w)u, w) = 0, so it follows that R{u, w)u = 0. A similar argument 
shows that R{w,u)w = 0. □ 

Proof of Theorem \1.4\ Since the manifold M is two dimensional, v = 
au + hw for some constants a and h. It follows that (15. ip becomes 

= {{Vu,R)iw, u)v, u) = b {{'VwR){ui, u)w, u) . 



However, by the proof of Theorem 15. 2[ {{W yjR){w,u)w,u) = 0. 
Therefore, fl5.ip is satisfied automatically. 

Since the Riemannian curvature R satisfies R{u, w)u = R{w, u)w = 
and M is 2-dimensional, the term G is simplified to 

G{u, v,w) = ^ {{VlR){v, u)v, u) 

+ \ {{V,V^R){w,u)v,u) + {{VlR){w,u)w,u) 
5 10 

= — {{VlR){w,u)w,u) + — {{V^VuR){w,u)w,u) 

ah c? 
+ y ((V„V^i?)(w, u)w, m) + ^ (iylR){w, u)w, u) . 

Let us extend u and v to vector fields V and VT, respectively, which 
are constant in a geodesic normal coordinate neighborhood. Since the 
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Gauss curvature at x vanishes and the covariant derivatives satisfy 



0, it follows that 



{{Vuyu,R)iw,u)w,u) = VuVw {R{W,U)W,U) 



= Vw^u{R{W,U)W,U) = {{V^VuR){w,u)w,u). 

X 

Therefore, the formula for G simplifies to 



(5.2) 



{{WlR){w,u)w,u) 



G{u, V, w 

+ {{V^VuR){w, u)w, u) + ^ {{VlR){w, u)w, u) 
5 lU 



Since the manifold M has non-negative Gauss curvature and has zero 
Gauss curvature at x. We have 



{{VlR){w, u)w, u) = Vl, {R{W, U)W, U) 



> 0. 



Therefore, if the quadratic in (15.21) satisfies G{u,v,w) > for all v, 
then the discriminant is non-positive and it follows that 

3 {{V^VuR){w, u)w, uf < 2 {{VlR){w, u)w, u) ((V^-R)(w, u)w, u) . 

□ 



6. Example 

In this section, we discuss the proof of Theorem 11.51 Recall that we 
consider the following Riemannian metric (-, ■) on M^: 

(6.1) (u, = e^^^^^u ■ t;. 

Let us denote the gradient and the Laplacian of the usual Euclidean 
metric by V and A, respectively. Let V and R be, respectively, the 
Levi-Civita connection and the Riemannian curvature of the Riemann- 
ian metric (■,■). 

Lemma 6.1. The Levi-Civita connection V and the Riemannian cur- 
vature R of the metric (■, ■) are given by 

'^d:,dy = Vdydx = fxdy + fydx, V d^O^ = —'^dy^y = fx^x " f yQy 

K = -(A/) ■ e-^f 

where K denotes the Gauss curvature with respect to the Riemannian 
metric (■, ■). 
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Proof. It follows immediately from, for instance, [21 Theorem 1.159]. 

□ 

Proposition 6.2. Let f{x, y) = x^y + ax^y"^ + xy^ . Then the Riemann- 



ian metric defined by ^6. 1\) has non-negative Gauss curvature if 

(6.2) a < -3. 

Proof. A computation shows that 

A/(a;, y) = 2ax^ + 12xy + 2ay'^. 

It follows from Lemma [6. II that the Gauss curvature is non- negative 
if and only if the quadratic 2ax^ + 12xy + 2ay'^ is non-positive. This, 
in turn, is equivalent to a < and 144 — 16a^ < 0. □ 

Proof of Theorem \1.5[ Assume that the cost c = d"^ satisfies the weak 
MTW conditon. By Proposition 16.21 and Theorem ll.il we have a < —3. 
Let U and W be two constant vector fields which are orthonormal 

= u and 



with respect to the Euclidean metric. Assume that U 



W 



(0,0) 



(0,0) 



w. Then 



(6.3) 



{R{U,W)U,W) = e^fR = 
= -e^^(A/) 
= -e^^ {2ax'^ + 12xy + 2ay'^) 
It follows from Lemma 16.11 that 



(0,0) 



Therefore, 

(6.4) (vlR{u,w)u,w 



(0,0) 



(0,0) 



Vl{R{U,W)U,W) 



(0,0) 



(0,0) 



0. 



Since Vg = e "^-^Vg for each smooth function g, it follows from 
Lemma 16.11 again that 



VuVg,uj = {VuVg,u). 
It follows from this, ([63D, and ([62]) that 

V^-R(m, w)u, = —4:{aul + 6uiU2 + aul). 
Similar calculations show that 

V'f^R{u, w)u, w) = —A{awl + 6wiW2 + aw^) 
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and 

V^Vu-R('U, w)u, w) = —A{auiWi + 3uiW2 + 3wiU2 + au2W2)- 



Since u and w are orthogonal, we can assume that w is given by 
w = —U2dx + uidy. It follows that 

3 {iyyjVuR){w, u)w, uf 

(6.5) - 2 {{VlR){w, u)w, u) {{VlR){w, u)w, u) 

= 16[(27 - 2a^)u^2 + (18 - 4:a^)ulul + (27 - 2a^)u^]. 

If 27 — 2a^ > 0, then we can set ui = and see that the above 
expression is positive for some u. This contradicts with Theorem 15.11 

Therefore, we have a < —\/^- Finally, we remark that the expression 



in (16. 5p is non-positive if a < D 



7. Proof of Theorem 15.11 



In this section, we will give the proof of Theorem 15.11 Let us first 
recall the notation that we are using. Let r i— f/s,t(''") be the parallel 
translation of the vector u along the geodesic r h-)- 7s,t(T) := exp(r(tf + 
sw)). Let T I— !■ Js,tij) be the Jacobi field defined along the geodesic -^s^t 
which satisfies the conditions Js,t(0) = m, Js,t{X) = 0; Js,t{,T) 7^ 
for all T in the interval (0, 1). 

First, it follows from (1) of Lemma [8.11 and (14. ip that 



dtdl{R{drl,J)d^^,U) 



s=t=0 



= {R{d%i, J)dtdrj, U) + {R{dtd^i, J)d%i, U) 
+ 2 {{dsR){dsdri. J)dtdri, U) + 2 ((9,i?) (9^9.7, J)dsdrl, U) 
+ 2 {R{dsdr^, dsJ)dtdri, U) + 2dt{{R{dAl. J)dsdrl, U)) 
+ 2 {R{dsdri, J)dtdri, dsU) + 2 (i?(9t9,7, J)9.9,7, t^) 

+ 2(i?(aia.7, J)aA7,5st/) 
By (2') of Lemma [8. 11 (1) of Lemma [8.21 Lemma [8. 3[ (1) of Lemma 



s=t=0 



9.H and (1) of Lemma [9. 3[ the above equation simplifies to 

dtdl{R{dr^,j)dr^,u) 



s=t=0 



(7.1) = 2r(l - r) ((V^i?)(w, u)v, u) + 2r(l - r) {{V^R){v, u)w, u) 
+ 2dt{{R{d,dr'y,J)dAl,U)) 



s=t=0 
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By (1) of Lemma (8 .21 (2) of Lemma [9.11 and (1) of Lemma [9.31 we 
also have 



dt{{R{dsd^^, J)d,dr7, U)) = {idtR){dsdr7, J)dsdr7, U) 



s=t=0 



s=t=0 



Therefore, it follows from (2') of Lemma [8.11 and Lemma [8.31 that 



dt{{R{d,drj, J)dA^, U)) = r(l - r) ((V,i?)(w, u)w, u) 



=t=o 



If we combine this with fl7.ip . then we have 



dtdl{R{drlJ)d^l.U) 



=t=0 



= 2r(l - r) {{Vu,R){w,u)v,u) + 2r(l - r) {{V u,R){y,u)w,u) 
+ 2t{1-t){{V^R){w,u)w,u). 

Therefore, by Theorem 14.11 we have 



dtMT^{u,tv,w 
1 



i=0 



4 



((V^-R)(u7, u)v, u) + ((V^-R)(t;, u)w, u) 



+ {{V^R){w,u)w,u) 



Finally, by taking covariant derivative of the property 



{R{w,u)v,u) = {R{v,u)w,u) 



of the Riemannian curvature R, we have 



dtWTW{u,tv,w 
1 



i=0 



2 ((V^i?)(w, u)v, u) + - {{VvR)iw, u)w, u) 
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By (1) of Lemma 18. 1[ (1) of Lemma 18.21 (1) and (2) of Lemma 19.11 
(1) of Lemma 19751 if we differentiate (14. ip twice with respect to t, then 



t=s=0 

= 2 {{d^^RMd^^, md^l, U) + 2 {R{dtdldri. J)dtd^l, U) 
+ 2 {R{dtdri, dU)dtdri, U) + 2 {R{dtdr^, J)dtd%^, U) 
+ 2 {R{dtdr-f, J)dtdri, dlU) + 4 {{dtdsR){dAi, J)dtdri, U) 
+ 4 {{dsR){dtd^^, dsJ)dtdrj, U)+4 {{mR){dAj, J)dA^, U) 
+ 4 {{dsR){dtdr-t, J)dtdr7, dsU) + 4 (/?(9,9,7, 9*5,7)9*9,7, U) 
+ 2 ((S^i?) (9,9,7, ^)9.9,7, U) + 2 (/2(929,9,7, J)9,9,7, U) 
+ 2 (i?(9,9,7, d^J)dsdri, U) + 2 (i?(9,9,7, J) 9*^9,9,7, U) 
+ 4 (i?(9,9,7, J)9*9,7, dtdsU) + 4 (i?(9i9,7, 9*9,7)9,9,7, U) 

+ 4 (i?(9t9,7, 9, J)9*9,7, 9,[/) + 4 (i?(9i9,7, J)9,9,7, dtd.U) 



By (2') of Lemma 18.11 Lemma 18.31 (3) of Lemma 19.11 Lemma 19.21 
(1) and (2)of Lemma the above equation simphfies to 



9,^9,^ (i?(9,7, J)drl.U) 



t=s=0 



= 2(1 - ry {{VlR){v, u)v, u) + 2(1 - ry {R{R{v, w)w, u)v, u) 

2t{t — l)(r — 2) 
H ^ {{R{v, R{w, u)w)v, u) + {R{w, R{y, u)v)w, u)) 

+ 4(1 - ry {{V,V^R)iw, u)v, u) + 4(1 - ry ((V,V^/?)(t;, u)w, u) 
+ 2(1 - ry {(ylR){w, u)w, u) + 2(1 - ry {R{R{v, w)v, u)w, u) 
+ 2(1 - ry {R{w, u)R{v, w)v, u) + 2(1 - ry {R{v, u)R{v , w)w, u) 
+ 2(1 - ry {R{w, u)v, R{v, w)u) + 2(1 - r)r^ {R{v, u)w, R{v, w)u) 
Ar{T- l)(r- 2) 



+ 



4r(r- l)(r + 1) 



((i?(w, R{w, u)v)w, u) + {R{w, R{w, u)v)v, u)) 
{{R{v, R{v, w)u)w, u) + {R{w, R{v, w)u)v, u)) 
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We simplify the above equation further by using the property of the 
Riemannian curvature R. 



t=s=0 



d^dUR{drl,J)dr7,U) 

= 2(1 - r)T^ {{\/lR){v,u)v,u) -4:{1 - t)t'^ {R{v,u)u, R{v,w)w) 

4r(r- l)(r-2) , 
+ ^ {R{v, u)v, R(w, u)w) 

+ 8(1 - ry ((V,V^i?)(w, u)v, u) 

+ 2(1 - r)r2 {{VlR){w, u)w, u) - 4(1 - {R{w, u)u, R{v, w)v) 

+ 2(1 - {R{w, u)v, R{v, w)u) + 2(1 - r)r^ {R{v, u)w, R{v, w)u) 

4r(T — 1)(t — 2) 

H ^ {{R{w, u)v, R{w, u)v) + {R{v, u)w, R{w, u)v)) 

3 

_ l)(r + l) 

^ -^^ ({R{w, u)v, R{v, w)u) + {R{v, u)w, R{v, w)u)) . 

If we integrate the above twice with respect to r and multiply by 
3/2, then we have 

d^MTW(u,tv,w) 

t=o 

1 1 
= Yo ~ 5 '^''^^ 

4 2 

+ — {R{v, u)v, R{w, u)w) + - ((V^V^i?)(w, u)v, u) 

-Lo O 

+ ^ {{^lR){w, u)w, u)-^ {R{w, u)u, R{v, w)v) 
4 

+ — {{R{w, u)v, R{w, u)v) + {R{v^ u)w, R{w, u)v)) 
io 

H — {{R{w, u)v, R{v, w)u) + {R{v, u)w, R{v, w)u)) . 
3 



8. Appendix 1: Lemmas for Natural Mechanical Actions 

In the two appendices, we give the proof of various lemmas used in 
the previous sections. The first appendix is devoted to those lemmas 
which are related to the natural mechanical actions. The rest of the 
lemmas needed only in the Riemannian case are done in the second 
appendix. Let us first recall our notations. Let u, v, and w be tangent 
vectors based at a point x which is a critical point of the potential V. 
Let T 7s,t(T) be the curve of least action with initial velocity tv + sw. 
Let T i-> Usfi{T) be the parallel translation of the vector u along the 
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curve r 7s,t(T). Let r Js,t{^) tic the Jacobi field defined along the 
curve T !->■ 7s,t(T) by the Jacobi equation 

dp + R{d + T-f, J)dr-i + Hessy ( J) = 

and satisfies the conditions Js,t{^) = Js,t{^) = 0, and Js,t{^) for 
all T in the interval (0, 1). Let r i->- v{t) be the solution to the initial 
value problem 

dp — —Y{.essVx{v) , V 



T=0 



0, drV 



V. 



T=0 



Similarly, let w be the solution to the above initial value problem with 
V replace by w. Let u be the solution to the boundary value problem 



dp = —iiessVx{u), u 



T=0 



M, u 



r=l 



0, u^Q if < r < 1. 



Lemma 8.1. The family of curves 7 satisfies the following: 

(1) 9,7 = 0, 

s=t=Q 

(2) % = V, ds-f 



s=t=0 



s=t=0 



W, 



In particular, if we are in the Riemannian case where the potential 
V = 0, then we have 



(2') 



5*7 



= TV, 



s=t=0 



TW. 



s=t=0 



Proof. Recall that 7 satisfies the Newton's equation dp — — VV^ with 



initial condition drj 



T=0 



tv + sw. Since x is a critical point of 



the potential V, it follows that 7 



= X is the solution to the 



s=t=0 



above initial value problem with s — t — 0. Therefore, (1) follows 
immediately from this. 

If we differentiate the Newton's equation with respect to t, then we 
have 

Rida, a,7)9r7 + dpa = dtdP = -}lessV^{dt-f). 
If we set s = t = and apply (1), wc have 

dptl . _ = -HessKr(% ). 



s=t=0 

We also have the initial conditions 

dtl 

t = S=T=0 

It follows that V — dt'f 

s=t=0 



= -Hess (9^7 
ions 

0, drdt^ 



s=t=0 



t=S=T=0 



□ 



Lemma 8.2. The family U of parallel vector fields satisfies 
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(1) d,U 

(2) dlU 

(3) d^,U 



dfU 



s=t=0 
s=t=0 
s=t=0 



s=t=0 



0, 



In particular, d^U 
where V = 0. 



JJ" R{drW, w)udT, 
R{drV, v)udT. 



s=t=0 



in the Riemannian case 



s=t=0 



Proof. Since the family is parallel, we have drU = 0. Therefore, we 
have 

= dsdrU = i?(9,7, dri)U + drdsU. 



If we evaluate at s = t = and apply (1) of Lemma [8.11 then we have 



drdsU 



0. 



Therefore, dgU 



s=t=0 

is constant in r and we have 



s=t=0 

dJJ 



s=t=0 



T = S=t=0 



dgU = 0. 



For the proof of (2), we apply drU = again and get 



d^dfU 



s=t=0 



dtdrdtU + R{drl,da)dtU 



s=t=0 



d^drU + dt{R{dr-f, da)u) + R{dri, da)dtU 



s=t=0 



By (1) and (2) of Lemma [8.11 the above equation becomes 



drd^U 



s=t=0 



d^drU + R{dtdrj,da)U 



s=t=0 



= R{drV, v)u. 
If we integrate with respect to the r-variable, then 



S=t=0,T=T 



+ / R{dT-v,v)udT 

T=S=t=0 Jq 



R{drV, v)udr. 



By (2') of Lemma I8.H v = tv in the Riemannian case. It follows 
from the skew symmetry of the Riemannian curvature that 



0. 



=i=0 



□ 
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Lemma 8.3. The family of Jacobi fields J satisfies 



J 



u 



s=t=0 



In particular, if we are in the Riemannian case where the potential 
V = 0, then we have 



J 



:i-t)u. 



=i=0 



Proof. Recall that the Jacobi equation is 

d'^J + R{dri, J)dri + Hess\/^( J) = 0. 

If we set s = t = and apply (1) of Lemma IHTTI then the above equation 
becomes 



(9^J + HessV;(J) 



0. 



s=t=0 



We also have the boundary conditions J 



u and J 



follows that J 



r=0 



0. It 



r=l 



U. 



s=t=0 



In the Riemannian case, V = and the above equation becomes 
0. If we combine this with the boundary conditions, we 

= (1-t)u. □ 



s=t=0 

have J 



s=t=0 

9. Appendix 2: Lemmas for the Riemannian case 

Let us first recall and specialize our notations used in the previous 
appendix to the Riemannian case. Let u, v, and w be tangent vectors 
at a point x and let r Js,t{T) := exp(r(tu + sw)) be the geodesic with 
initial velocity tv + sw. Let r t— )• Us,tiT) be the parallel translation of 
the vector u along the curve r i— )■ 7s,f(r). Let r i-> Js,t{j) be the Jacobi 
field defined along the curve r i-> 7s,f(T) which satisfies the conditions 
Js,tiS^) = Js^iX) = 0, and Js,tij) 7^ for all r in the interval (0, 1). 

Lemma 9.1. The family of geodesies 7 satisfies the following: 



(1) dld^^ 

(2) dtOAi 

(3) dtd%^ 

(4) 929,9.7 



s=t=0 
s=t=Q 
t=0 



s=t=0 

= dsdtdr-f 



0, 



0. 



s=t=0 



t'^R{v,w)w, dgd^dr'y 



s=t=0 
s=t=0 



= t"^ R{v , w)v , d'^dtdr'j 

s=t=0 

Proof. For (1), we have 

(9.1) 9,29.7 = ^r^h + Rids^, 5r7)<9.7- 



= T^R{w, v)v, 
T^R{w, v)w. 
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Since 7 



i=0 



exp{sTw) is a geodesic in the variable s, we have 



0. 



=i=0 



Therefore, if we set s = t = 0, then (19. ip becomes 



=i=0 



=t=0 



Finally, if we apply (1) of Lemma then we obtain (1). 
For (2), we have 



drdtdsdr'y 

s=t=0 



s=t=Q 



s=t=0 



If we apply (1) of Lemma [8.11 then the above equation becomes 



drdtdsdri 



=t=0 



s=t=0 



Since 7 is a geodesic for each t and s (i.e. 9^7 = 0), it follows that 



drdtdsdr'y 



0. 



Therefore, dtdsdr'y 
dtdsdrl 



s=t=0 

is constant in r and we have 



s=t=0 

dtdsdri 



s=t=0 

For (3), we have 

drdtdldri 



=t=o 



dtds{tv + sw) 



0. 



s=t=0 



s=t=0 



s=t=0 



By (1) of Lemma [8. H the above equation becomes 



drdtdldr-i 



s=t=0 



dtdrdldr^i 



=i=0 



If we apply (1) of Lemma [8?T] again to the above equation, then we 
obtain 



drdtdldr-i 



s=t=Q 



dtdsdrdsdri + R{dtdri, dsi)dsdri 



s=t=0 



s=t=0 
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By (2') of Lemma [8.11 we have 
drdtd'^^dr-f 



s=t=0 



dtdgdrdgdr^ + tR{v, w)w 



s=t=0 



Since r h-)- 7 is a geodesic, the above equation becomes 



s=t=0 



dtds{R{dri, ds'j)d-r'j) + tR{v, w)w 



s=t=0 



If we apply (1) and (2') of Lemma [8.11 then we have 



s=t=0 



= R{dtdr'y, dsi)dsdri + R{dsdr'y, dsi)dtdri + tR{v, w)w 
= 2tR{v, w)w. 
If we integrate the above equation in r, then we get 



s=t=0 



dtdldri 



s=t=0 



t'^R{v, w)w + dtdldr'^ 



T=S = t=0 



t'^R{v, w)w + dtdg{tv + sw] 
t'^R{v, w)w. 



T=S = t=0 



This finishes the proof of (3). 

For (4), we first apply (1) of Lemma El 



drd^dsdr'y 



s=t=0 



R{dr'y, dt'y)dtdsdr'j + dtdrdtdsdr'j 



s=t=0 



dtdrdtdsdri 



s=t=0 



By (1) of Lemma [8.11 again, the above equation becomes 



drd^dsdr-f 



=t=o 



s=t=0 



=t=0 



By (2') of Lemma [8. 11 we have 



s=t=0 



d'^drdgdr'y + tR{v, v)w 



=f=0 



d'fdrdsdr'y 



=t=0 



Since r h-)- 7 is a geodesic, we get 



drdfdsdr-f 



s=t=0 



s=t=0 



dKR{drl,dsj)drl) 



s=t=0 



COMPUTABLE NECESSARY CONDITIONS FOR SMOOTH OPTIMAL MAP 23 

By (1) and (2') of Lemma [8 -H the above equation becomes 



s=t=0 



s=t=0 



2tR{v, w)v 



s=t=0 



Finally, if we integrate the above equation in r, then we obtain 



d^dsdr-i 



s=t=0 



T = S=t=0 



t'^R{v, w)v + dldsdr'^ 

t'^R{v, w)v + dtds{tv + sw) 
r'^R{v, w)v. 



s=t=0 



Lemma 9.2. The family of parallel vector fields U satisfies 



T 



R{v,w)u, dgdtU 



s=t=o 2 
Proof. By (1) of Lemma [8. we have 



T 



s=t=0 



R{w, v)u. 



drdtdsU 



s=t=0 



dtdrdsU + R{dri, da)dsU 
dtdrdsU 



s=t=0 



s=t=0 



If we apply (1) of Lemma [8.11 again, then the above becomes 



drdtdsU = mdrU + dt{R{drl, dsl)U) 

s=t=0 

= dtdsdrU + R{dtdri,dsi)U 
Since r i— )• f/ is a parallel vector field, we get 



s=t=0 



s=t=0 



drdtdsU 



s=t=0 



R{dtd^j,d,j)U 



s=t=0 



If we apply (1) and (2') of Lemma [8 .11 then we have 



drdtdsU 



s=t=0 



tR{v, w)u. 



Since U 



□ 



r=0 



u, we can integrate the above equation in r and obtain 



T 



s=t=0 



R{v,w)u + dtdsU 



T=S = t=0 



-R{v, w)u. 



□ 



Lemma 9.3. The family of Jacobi fields J satisfies 



s=t=0 



dtJ 



s=t=0 



0, 
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T(T-l)(r-2) 



r(r-l)(r-2) 



, u)v, 
R{w, u)w, 



=i=0 



Proof. The family J satisfies the Jacobi equation 

dlJ + R{dri,J)dri = 



and the boundary conditions J 



u and J 



T=0 



0. 



r=l 



If we differentiate the Jacobi equation with respect to s and apply 
(1) of Lemma [8. H then we have 

dsdlj = -ds{R{drl,J)drl) 

It follows that 
= d.dlJ 

s=t=0 



drdsdrJ + R{ds'y, dr^)drJ 



s=t=Q 



s=t=0 



By (1) of Lemma [8 -H we have 



0. 



s=t=0 



T=0 



dsJ 







T=l 



This together with the boundary conditions dgJ 
give dgJ =0. This finishes the proof of (1). 

s=t=0 

For (2), we differentiate the Jacobi equation with respect to t twice 
and apply Lemma [SHI 



dplJ 



s=t=0 



-d^^{R{drl,J)drl) 
-2R{dtdrl, J)dtdrl 



s=t=0 



s=t=0 



Therefore, by (2') of Lemma [8.11 and Lemma [8.3[ we have 



(9.2) 



s=t=0 



-2{1 -t)R{v,u)v. 



On the other hand, by (1) of Lemma [8.11 we have 



dfdp 



s=t=0 



dtdrdtdrJ + dt{R{da, dri)drJ) 
dtdrdtdrJ + R{da, dtdri)drJ 



s=t=0 



s=t=0 
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By (2') of Lemma [8.11 the above equation becomes 

dtdrdtdrJ + tR{v, v)drJ 



=t=0 



s=t=0 



dtdrdtdrJ 



s=t=0 



drd^drJ + -R(<9i7, dr'^)dtdrJ 



s=t=0 



If we apply (1) of Lemma [8.11 then we have 



s=t=0 



drdfdrJ 



s=t=0 



drdtdrdtJ + drdt{R{da, dri)J) 
drdtdrdtJ 



s=t=0 

By (2') of Lemma [8. 11 we have 



s=t=0 

+ dr{R{da,dtdrl)J) 



s=t=0 



s=t=0 



drdtdrdtJ 
drdtdrdtJ 



s=t=0 



s=t=0 



+ dr{R{TV,v)J) 



s=t=0 



By applying (1) of Lemma [8.11 the above equation becomes 



d^J 



=t=o 



d',dtJ + dr{R{dt^,dri)dtJ) 



s=t=0 



dldlJ 



s=t=0 



T=0 



du 



T=l 



If we combine this with fl9.2p . then we have 

d'^d'^ J = -2{1 - t)R{v,u)v. 

This together with the boundary conditions dfJ 
gives us 

This finishes the proof of (2). The proof of (3) is the same with v and 
w interchange. 

For (4), we apply (1) of Lemma [8.11 and get 

d%d,J 

s=t=0 

= drdtdrdsJ + dr{R{dri, dti)dsJ) 



S=t=Q 



drdtdrdsJ 



s=t=0 
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If we apply again (1) of Lemma ISTTl we get 

dldtdj 

s=t=0 



s=t=0 



dtd%J + R{dr7,da)drdsJ 

dtd%J 

s=t=0 

dtdrdsdrJ + dtdr{R{dr^, dsi)J) 



s=t=0 



Using the fact that r i— t- 7 is a geodesic and applying (1) of Lemma 
8.11 we get 



dldfAJ 



s=t=0 



= dtdrdsdrJ + R{dtdri, drdsi)J) + R{dtdri, dsi)drJ 
If we apply (2') of Lemma [8. II and Lemma [8.31 then we have 
dldtdsJ 

s=t=0 

+ {1 -2t)R{v,w)u 



t=0 



s=t=0 

dtdrdsdrJ 



s=t=0 



= dtdsdlJ + dt{R{dr^, dsi)drJ) + (1 - 2t)R{v, w)u. 

s=t=0 

By (1), (2') of Lemma [8. 11 and Lemma [8 .31 we get 



d^rdtdsJ 



s=t=0 



dtdsdlJ 



s=t=0 



+ (1 - 3t)R{v,w)u. 



If we apply the Jacobi equation and (1) of Lemma [8. ![ then 
d^dtdsJ 

s=t=0 

= -dtds{R{dri, J)dri) + (1 - 3r)i?(t;, w)u 

s=t=0 



s=t=0 



= -R{dsdr-f, J)dtdrj) - R{dtdr-f, J)dsdrj) 
+ (1 - 3t)R{v,w)u 
By (2') of Lemma [8. II and Lemma [8.31 we have 

dldtdsJ 

s=t=0 

= -(1 - t){R{w, u)v + R{v, u)w) + (1 - 3t)R{v, w)u. 
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Finally, by the first Bianchi identity, we have 

dldtdsJ = -2(1 - t)R{w, u)v - 2tR{v, w)u. 

s=t=0 



= dtdsJ 

T=0 



T=l 



0, we 



Using this and the boundary conditions dtdgJ 
have 

dtdsJ = ~ [(r - 2)R{w,u)v - (r + l)R{v, w)u]. 

□ 
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